Anomalous stabilization in a spin-transfer system at high spin polarization 
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Switching diagrams of nanoscale ferromagnets driven by a spin-transfer torque are studied in the 
macrospin approximation. We consider a disk-shaped free layer with in-plane easy axis and external 
magnetic field directed in-plane at 90° to that axis. It is shown that this configuration is sensitive 
to the angular dependence of the spin-transfer efficiency factor and can be used to experimentally 
distinguish between different forms of g{0), in particular between the original Slonczewski form and 
the constant g approximation. The difference in switching diagrams is especially pronounced at 
large spin polarizations, with the Slonczewski case exhibiting an anomalous region. 



I. INTRODUCTION 

Spin polarized electric currents have been successfully 
used to switch the magnetization direction of nanoscale 
ferromagnetic layers via the spin transfer effeciji'^'^i^i^i^. 
One of the questions of current-induced dynamics is the 
dependence of spin-transfer efficiency, or Slonczewski fac- 
tor g, on the angle between the polarization of incoming 
spin current and the magnetization directionii^iiSiiiS. Such 
a dependence can be essential, and for example leads to 
the asymmetry between the positive and negative switch- 
ing currents. However, there is still a lack of experimental 
tests for the precise functional form of efficiency factor. 
It is expected that angular dependence of g will become 
more important at high spin polarizations where the con- 
stant efficiency approximation can fail, while constant g 
can still be in good agreement with experimental results 
at low spin polarizationiiii^ii^iii^. 

Here we perform stability analysis for the equilibrium 
configurations of a bilayer spin-transfer device using the 
Slonczewski form for the efficiency factor and compare it 
with a similar analysis that uses the constant efficiency 
approximation. We observe that the switching diagram 
for the Slonczewski case displays a stability region and 
precessional states that are absent in the constant effi- 
ciency case. These anomalous regions become larger as 
the spin polarization increases. Our results may moti- 
vate further experimental efforts to directly measure the 
functional form of the efficiency factor at high spin po- 
larizations. 



II. MACROSPIN DESCRIPTION OF THE 
DEVICE 

A typical device used to study spin-transfer effect is 
a nanopillar, with two layers of ferromagnetic material 
separated by a normal paramagnetic metal (see fig [TJa). 
The magnetization of one layer (polarizer) is fixed and 
oriented along a unit vector s, while the magnetization of 
the other (free layer), M = Mn, rotates and is described 
in the macrospin approximation by the Landau-Lifshitz- 
Gilbert (LLG) equation including the Slonczewski spin 
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g{9,P)[nx (sxn)] + a[nxri], 

(1) 

where 7 is the gyromagnetic ratio, E{n) is the magnetic 
energy of the free layer, and a is the Gilbert damping 
constant. The strength of the spin-torque is character- 
ized by the efficiency factor, g{d,P), which depends on 
the angle 9 between the magnetizations of the polarizer 
and the free layer, and the degree of current spin po- 
larization P £ [0,1]. In general the functional form of 
g{9, P) is material and geometry-dependent^i^'^'i^. Here 
we will compare Slonczewski's^ form 



9{0,P) = 
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(1 + Pf/AP^/^ and the 



with = 3 - 4//p, /p 
g{9, P)=const approximation. 

The magnetic energy of the free layer includes con- 
tributions from the intrinsic anisotropy (easy axis 
anisotropy with strength Ha and direction a), shape 
anisotropy (easy plane anisotropy with normal vector 
p) and the interaction energy with an external magnetic 
field H. Equation ([1]) can be written as 

n — T(n) + cm x 'r(n) , (3) 

where we have rescaled the time as T = i/(l + a^), and 
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FIG. 1: a) Typical nanopillar device with free layer on the 
top and polarizer at the bottom, b) In-plane magnetic field 
configuration. Two stable directions of the free layer magne- 
tization are labeled as N and S. 
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T is defined as 

T(n) = — V£(n) X n + L0ig{6, P) n x (s x n), 



e(n) 



-(p.n)2-^(a.n)2-c.H(h-n) 



(4) 



The newly defined constants are related to the already 
introduced parameters according to 



UJH = iH, 



(5) 
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All of them have dimensions of frequency making the 
comparison between the terms of different origin straight- 
forward. In accord with experimental situations it is as- 
sumed that uji <^ UJp. 

We study a device with an in-plane easy axis and the 
in-plane magnetic field perpendicular to it (see fig. [lib). 
Choosing the system of coordinates § = a = Gz, h = Gx, 
p = Gy, we obtain, from equation ([¥]), the components of 
T in spherical coordinates 



T0 = i sin 29{uLip s\v? -I- Wa) - luh 



cos a cos ( 



Tg = ^ sin 9 sin 20 — uoh sin — uiig{cos 0) sin 9. 



(6) 



The equilibrium directions of the magnetization n corre- 
spond to the solutions of the equation T(n) = 0. Here we 
consider the two in-plane equilibrium points. At u)i = 0, 
ijJh = these are the north (N) and the south (S) poles. 
For u)j — 0, LdH they shift and approach the direc- 
tion of magnetic field, finally merging at ojh = Wo. The 
shifted equilibrium points are still labeled by N and S 
(Fig.Hh). 

Stability of an equilibrium can be checked by expand- 
ing T in angular deviations 59, S(f>, and writing equation 
([3]) in an approximate form 
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The equilibrium is stable when the real parts of both 
eigenvalues of D are negative, or equivalently when ma- 
trix D satisfies TrD < and detD > at the equilibrium. 



III. STABILITY REGIONS 

The modified positions of the N and S equilibria for 
LOH ^ 0, ^ are given by 



sm^AT s = h O — 
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(8) 



for Q < ujh < uja (with gN,s = 9iON,s))- The angle 
9 and the magnetic field strength ujh have a one-to-one 
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FIG. 2: Switching diagram for the g =const approximation. 
The value of g is chosen as the average value of Slonczewski's 
g{6) used in Fig. [S] Other parameters are set to uia/ujp = 
0.01, P = 0.7, and a = 0.01. Stability regions for the N 
and S equilibria (see text) overlap forming the bistable region 
marked as B. 



correspondence and can be used interchangeably. The 
trace of D-matrix at these points can be found as 



TrD = -uJi 



g{9)cos9+-{g{9)sml 



(9) 



Approximations 
form 



a[ujp cos 20 -I- {ujp sin^ + u;a){l + cos^ 9)] 

give the stability condition in the 
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The determinant 
1 



1 



detD = uja{ujp + uja) cos-' 9 + 0{ujf), (11) 



in the small current regime remains positive, so it does 
not play any role in the stability analysis in this case. In 
contrast, the trace TrD is more sensitive and can change 
sign as the current is varied. Moreover, the explicit ap- 
pearance of g'{9) in the formula leads to important dif- 
ferences in the switching diagrams for different forms of 
g{9). In the case of constant g-factor the stability condi- 
tion for N- and S-equilibira can be written as 



uji ^ =pa 



UJp+U}a{2 - {UJH/Ujgf) 

2g^l - [ujH/uJa)^ 



(12) 



where >, — (<, +) corresponds to the N (S) stability 
region (See Fig. The switching current exhibits the 
1 / cos 9 divergence reported in the experiments for this 



regime^. 



For the Slonczewski g-factor, the condition of stability 
for the N point is 



3 



(Ol 
(Op 

0.05 






- 







-0.05 


1 






0.2 0.4 0.6 


0-8 COh 1 

co„ 



FIG. 3: Switching diagrams for the Slonczewski's form of g- 
factor. Other parameters are the same as in Fig. (2] Regions 
of stability for north (N) and south (S) poles overlapping in 
the bistable region (B), and the region of precessional states 
(P) are shown. The onset of the anomalous stability behavior 
occurs at a field uj'^/lop — 0.76. 



In the "anomalous" regime loh > ^ current of pos- 
itive polarity stabilizes both N and S points, while the 
application of a sufficiently large negative current desta- 
bilizes both points (see Fig.[21), moreover, there is a region 
in which none of the equilibria are stable, suggesting the 
existence of precessional motion. The value of lo*]^ be- 
comes smaller as the polarization increases. In the limit 
P — > 1 it becomes zero, so that the anomalous region 
fills all the switching diagram. In other words, as the 
polarization becomes larger the differences between the 
g-constant approximation and the Slonczewski form be- 
come quite dramatic. The position 9c of the S point at 
Lo*u is shown in Fig. ID 

Substantial difference between the switching diagrams 
at large spin polarizations found in this study underscores 
the necessity of developing new experiments capable of 
determining the g{9) dependence. It also suggests that 
in the regime of large spin polarization the behavior of 
spin-transfer devices may experience qualitative changes. 
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whereas for the S point the condition becomes 



, (13) 



25S\/1 - {uJH/^^aY - glfp{uJH/l^a) 
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where < (>) is the condition for lu < lo*^ (uj > Lo*fj), 
designates the field for which the denominator of equa- 
tion (|14p becomes zero, and determines the onset of an 
stability behavior completely absent in the g-constant 
case (Fig. [3]) . This field, or equivalently the angle char- 
acterizing the S point, depends only on polarization P 
and can be found from 
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FIG. 4: Critical angle 6c for the onset of the anomalous sta- 
bilization as function of the polarization 
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